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Abstract 

Coupled map lattices have been widely used as models in several fields 
of physics, such as chaotic strings, turbulence, and phase transitions, as 
well as in other disciplines, such as biology (ecology, evolution) and infor- 
mation processing. This paper investigates properties of periodic orbits 
in two coupled Tchebyscheff maps. Then zeta function cycle expansions 
are used to compute dynamical averages appearing in Beck's theory of 
chaotic strings. The results show close agreement with direct simulation 
for most values of the coupling parameter, and yield information about 
the system complementary to that of direct simulation. 

1 Introduction 

Coupled map lattices (CML) are dynamical systems vi/ith discrete space and 
time and continuous state variables. Their most common form in applications 
is the diffusively coupled model 

1 = (1 - + (a/2)[.9($(:+i') + gi^t'^)] (1) 

where n is a discrete time step, z is a lattice point with a periodic boundary 
condition, /($) the local map, g{^) the coupling and a a continuous parameter. 
They prove particularly useful in the investigation of chaotic dynamics of spa- 
tially extended systems (spatiotemporal chaos), but have also been employed 
in chemical, biological and engineering modelling, as well as pattern formation. 
The discrete structure in space and time can either be the approximation of 
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a less tractable continuous system (as in turbulence), or may arise naturally 
in applications (as in population dynamics of discrete colonies of organisms). 
Refs. m 121 El El give a good overview of early work to 1998 including the above 
applications; recent developments include advances in the mathematics of the 
relevant Banach spaces and transfer operators jSl jS] as well as ever widening 
apphcations; examples include ecology 0111121; neurons ^lEl, traffic flow ^2 
cryptography ^1 ^1 and field theory |15[ [TBI ■ 

From the point of view of understanding the dynamics of coupled chaotic 
maps, the most common CML studied are where the local dynamics / is quadratic 
(the "logistic" map); this class of one dimensional maps is now reasonably well 
understood. The simplest example of a quadratic map is the second degree 
Tchebyscheff polynomial f{^) — T'2(<I>) = 2<I>^ — 1 which is well known to be 
conjugated to a piecewise linear map with complete binary symbolic dynamics, 
and hence exactly solvable. 

In this paper we consider diffusively coupled Tchebyscheff maps (see Sec. 
2.2 for details) of degree 2 and 3. These have been previously used by Beck 
| 15l in some interesting but unconfirmed ideas concerning dynamical gener- 
ation of noise for stochastic quantisation leading to predictions of parameters 
in the standard model of particle physics. One of the authors ^71 subsequently 
showed that four of Beck's six CML models exhibit stable synchronised time pe- 
riodic states at arbitrarily small values of the coupling parameter a, albeit with 
small basins of attraction. The motivation for considering these models derives 
not only from possible applications, but also theoretical simplicity in that the 
uncoupled state is fully chaotic and exactly solvable, the small coupling limit has 
interesting structure, and the equations (being low degree polynomials) allow 
many exact calculations. 

Here we apply classical periodic orbit theory jl8[ [T^ to a lattice of two cou- 
pled Tchebyscheff maps. This theory allows very precise computations of long 
time properties (averages, Lyapunov exponents, dimensions, decay of correla- 
tions, etc.) of chaotic systems, but has been mostly restricted to low dimensional 
systems in the past. This paper contains the first steps towards developing a high 
dimensional periodic orbit theory, making use of the properties of Tchebyscheff 
maps to locate a large number of periodic orbits. If the Beck theory turns out 
to be true, the insight gained regarding these CML will be invaluable in devel- 
oping periodic orbit based computation methods which are potentially far more 
precise than direct simulation, and which are required in particle physics. More 
generally, a periodic orbit theory of CML will allow a far greater precision and 
insight into the many spatiotemporal systems discussed at the beginning than 
is currently available using direct methods. 

In order to achieve our immediate goal of computing Beck's observables on 
a two site lattice, we first need to look for periodic orbits of Q by using appro- 
priate numerical methods and study some of their properties, as described in 
the following sections. The locations, the stability and the number of periodic 
orbits of such maps depend on the value of a and current results do not allow 
prediction of these changes in detail. As we discuss in section 2, unexpected 
phenomena take place as soon as the coupling parameter a is increased from 
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0, for example a few periodic points suddenly move off the region G [—1,1], 
which we consider for our analysis. As a is increased from, periodic points drift 
on the domain G [—1, 1], moving into the complex plane and then reappear- 
ing on the real axis as a is increased further. The stability of periodic orbits 
is also heavily affected by the value of the coupling and stable and/or almost 
stable orbits appear for certain a's. Such behaviour has important consequences 
for the convergence of cycle expansions 1^, as discussed in section 3. In fact 
both the rate and the kind of convergence of any average we compute on the 
lattice strongly depends on the value of the coupling parameter. Still in the 
same section we address the applications of equations (0)1 to the theory of Beck. 
As pointed out by one of the authors J7] , the system does not show ergod- 
icity for all values of a, therefore direct simulation is not always guaranteed to 
be reliable in computing averages over the lattice. That is one reason why we 
investigate the efficiency of averaging over the periodic points. When it is possi- 
ble to achieve a significant precision for the zeta function (see sec. 3.1), we have 
computed the average of some string potentials defined in sec. 3.2, and then 
compared the results with direct simulation on the same lattice. Conclusions 
and outlook are presented in section 4. 

2 Periodic solutions 
2.1 Finding periodic orbits 

The first problem is to find periodic solutions to the map l@J up to the highest 
possible period. A number of numerical methods can be employed for that 
purpose. In our case the quickest one seems to be the inverse iteration method, 
as follows: Let x' = f{x) be a map whose cycles are all expanding^, then it 
follows that its inverse x' = f~^{x) has only stable periodic orbits, which can 
easily be found by iteration. A different symbol sequence (where symbols are 
associated with branches of the multivalued inverse) is to be used for every orbit 
|19| . Unfortunately this method cannot be applied a priori, unless we are certain 
that the map has only expanding cycles. It turns out that the system Q has got 
a few periodic orbits whose Jacobian has one stable and one unstable eigenvalue 
(see sec 2.3). This makes it impossible for inverse iteration to find them, since 
we cannot guess initial conditions that lie right on the unstable manifold. 

A different approach to the problem is that to find the zeros of the equation 
/"(x) — a; = 0, where n is the period in question. The multipoint Newton's 
method is believed to be one of the best algorithms for that purpose. It 
has the advantage of converging relatively fast to the solution, provided that 
sufficiently accurate initial conditions are given. A cycle of period n is a root of 

^they have no marginal or contracting directions 
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Figure 1: T2 and T3 polynomials as defined in ||2J). 



the vector function F: 





1 Xi 


] 




( X, 


- f{Xn) 


\ 


F{x) = F 


X2 






X2 


-fixi) 








) 






' f{Xn-l) 


J 



The iteration we need in order to approach the roots has the form 
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Therefore, after a certain number of iterations we expect to find one root. Of 
course the method has to be performed every time with different initial condi- 
tions, until we are sure all the solutions have been found. The issue concerning 
the choice of initial conditions is a delicate one, because in principle it is difficult 
to guess the initial points on a 2-dimensional plane accurately. Nevertheless, it 
is useful to remember that we are dealing with weakly coupled maps (a <C 1). 
For that reason we expect the periodic points to be relatively close to those of 
the uncoupled system (a = 0), that can be found analytically (see sec. 2.2) and 
hence exactly. 



2.2 Uncoupled system 

The spatially uncoupled dynamics of Tchebyscheff maps is given by 

$;+l=rAr(<i>;) 



(2) 
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where T/v is the Nth Tchebyscheff polynomial, defined in general as 

rjv($) = cos (N arccos $) (3) 

of which we consider the first two, namely 

T2($) = 2$2 _ 1 

T3($) = 4$3 -3$ 
The general (coupled) Tchebyscheff maps have the form 

= (1 - a)T^{^l,) + s^(T^($-i) + T'^{K+^)) (4) 

where a is the string coupling parameter, s can be ±1 and b can take on the 
values 1 and 0, determining respectively advanced coupling (T^($) — Tn{^), 
"A") or backward coupling (T^($) = 4>, "B"). We have then six models 2A+, 
2A~ , 2B~, 3A+ and 35"'". The other two = 3 maps are equivalent to 

these by a trivial symmetry. 
The model becomes 

Xn+i = (1 - a)TN{xn) + saT^{yn) (5) 
Un+i = (1 - a)TN{yn) + saT%{xn) (6) 

on a two site lattice. In this section we consider the uncoupled case a = 0, in 
which the system H5I6|I reduces to 

x'^Tm{x) v'^TNiy) (7) 

This is useful as there is a close correspondence between orbits at a = and 
orbits with a > 0. This correspondence is used to find orbits in the coupled 
dynamics using the uncoupled orbits as starting points for the multipoint New- 
ton algorithm. The correspondence is not exact — the A^^"p degree polynomial 
equation giving the periodic condition always has this number of (possibly de- 
generate) solutions, however for a > some of these solutions are not in the 
range [—1, 1] or are complex, and hence not relevant to the dynamics. This is 
discussed further in the next section. 

We begin first with some definitions. A periodic point of length Up is a 
solution of f"''{x) — x; a.s discussed above, there are at most A^^"?" of them. 
A periodic orbit is the orbit {f-'{x)},j = 0, 1, . . . Up — 1 of a periodic point. A 
prime orbit is a periodic orbit which is not the repeat of a shorter orbit, that is, 
(x) = X only if j is a multiple of rip. {x, y) is a periodic point of length rip of the 
uncoupled system if and only if x and y are both periodic points of length rip of 
the map Tm- It is well known JH| that the Tjv maps are conjugate to piecewise 
linear maps with complete N-avy symbolic dynamics, by the transformation 
X = cos tt9. This leads to exact expressions for all the periodic points of length 
rip, namely x = cos[2j7r/(A^"p ± 1)] with < j < A^"p/2. From these we 
enumerate all the prime orbits of the uncoupled two dimensional dynamics. 
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eliminating cyclic permutations which correspond to points on the same orbit, 
and repeats of shorter orbits. 

As an example we consider N — rip — 2. Let symbol (1) denote x € 
[—1,0] (respectively a; € (0, 1]); in general there are N symbols. In the symbolic 
representation of periodic points, the spatial variable {i in Eq. (1)) will increase 
to the right, and the temporal variable (n in Eq. (1)) will increase downwards. 
The periodic points of the one dimensional map are 



leading to 16 periodic points in the two dimensional map, 
orbits 

00 00 10 01 00 

01 ' 11 ' 11 ' 10 ' 10 ' 



their cyclic permutations 



11 

' 1 ' ' 1 



00 01 10 11 
00 ' 01 ' 10 ' 11 



classified as six prime 

01 
11 



11 
01 



01 11 11 10 10 
00 ' 00 ' 10 ' 01 ' 00 



and four repeated period one orbits 



2.3 Symmetry 

Until now we have not considered the symmetry of the system. This is partic- 
ularly important for cycle expansions |19l I2()| . for which a correct treatment of 
symmetry significantly accelerates convergence of the expansions. An uncoupled 
lattice has symmetry group Sx x where X is the number of lattice sites 
(equal to two in this paper), Sx is the symmetric group, of order X!, corre- 
sponding to permutations of the lattice points, and G is the symmetry group of 
the dynamics. A symmetrically coupled lattice with X > 2 has symmetry group 
-Dx, the dihedral group, of order 2X, corresponding to spatial cyclic permuta- 
tions and reflections of the lattice sites. In the case X = 2 the cyclic shift and 
reflection are equivalent, so the group is C2, the cyclic group. Remnants of the 
internal symmetry group G may or may not survive, depending on the nature 
of the coupling. For A" > 3 the coupling breaks the permutation symmetry, 
reducing the order of the symmetry group (ignoring internal symmetries) , from 
A! to 2X Elj. 

Each symmetry operation permutes periodic points. If the image of a peri- 
odic point lies in the same periodic orbit, that orbit is called symmetric, and is 
invariant under a cyclic subgroup of symmetry operations. If it does not, the 
new orbit is equivalent to the first under the symmetry, and shares many of the 
same properties, in particular stability and averages invariant under the sym- 
metry. It is also possible for a periodic point to be a fixed point of the symmetry 
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transformation, in which case its orbit is caUcd a boundary orbit. In the pres- 
ence of symmetry, the best strategy is desymmetrisation, in which all images of 
a phase space point are considered equivalent, and so may be described by a rep- 
resentative point on a subset of the original phase space called the fundamental 
domain. Equivalent orbits in the full dynamics correspond to a single orbit in 
the desymmetrised dynamics, and symmetric orbits in the full dynamics corre- 
spond to repeats of shorter orbits in the desymmetrised dynamics. Thus the 
enumeration of prime orbits up to a given period in the desymmetrised dynamics 
contains much more information than the enumeration up to the same period 
in the original dynamics, hence the accelerated convergence of cycle expansions. 

We continue with the previous example, considering the N = 2 maps, which 
have no internal symmetry. Thus the symmetry group is the cyclic group C2, 
with the nontrivial element corresponding to {x,y) {y,x). The fundamental 
domain can be taken as a; > y, from which the full domain is obtained by taking 
symmetry operations. The boundary satisfies x = y. Now the prime orbits of 
period two in the desymmetrised dynamics are classified as follows. There are 
two asymmetric orbits of period two in the original dynamics, 

00 10 
10 ' 11 

one boundary orbit of original period two 

00 
11 

and three symmetric orbits of original period four 

00 01 10 

01 01 11 
00 ' 10 ' 10 
10 10 11 

Note that the total number (six) is the same, but the information content is 
higher, since all the original period two orbits are included as period one or two, 
and also some orbits of original period four. 

In the case iV = 3, the symmetry group of the interacting case is C|, since 
there is an internal symmetry x —x, which in the coupled map survives 
as {x,y) — > {—x,—y). This is in addition to the permutation {x,y) — > {y,x). 
Thus the fundamental domain can be taken as x > \y\. Orbits may now be 
classified as asymmetric, symmetric under {x,y) {y,x), symmetric under 
{x,y) — > {—x,—y), symmetric under {x,y) — > {—y,—x), boundary with x = y, 
boundary with x = —y and (if it occurs) boundary with x = y = 0. 

2.4 Behaviour in the coupled case 

In this section we illustrate the behaviour of the map periodic solutions when 
the coupling parameter a is in the range [0,0.1]. In particular, we describe two 
phenomena, that is the way and stability and number of cycles change with a. 
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Stability is expressed by the real parts'^ Ai,A2 of the eigenvalues of the prod- 
uct of the Jacobian matrices 



J(x) 



( 




) 



of each point of one cycle, where f{x, y) and f{y, x) are the set of equation H5I6|I . 
A is said to be stable when |A| < 1. When a = 0, the systems is uncoupled 
and the maps have no stable cycles. Since the eigenvalues are continuous 
functions of a, we expect no stable eigenvalue to show up for very small a at 
fixed orbit period. As for stability, we analyse the trend of the maximum and 
minimum eigenvalues found for each period. In general, we consider the product 



so that in the 2-d case it is either A1A2 or just the unstable A. Stable orbits are 
not used for cycle expansions, hence they are not considered for this analysis. 
The graphs in Figure |21 show that both the maximum and the minimum eigen- 
values increase hyperbolically (i.e. exponentially) with respect to the period rip 
of the cycles when a — Q and still up to a ~ 10^^, but the trend of the minimum 
A's changes as long as a grows, so that their increase gets slower and slower. 
Besides, as we will explain later in this section, the few cycles with |A| = 4^"'p 
move off the fundamental domain when a 7^ 0, so that the maximum eigenvalues 
take the form |A| < 4"''.^ When a — 0.1, the minimum eigenvalues are per- 
haps best fitted by a power law, rather than an exponential. Such a behaviour 
is common to all the coupled T2 and T3 maps. Figure [3 also shows that the 
trend of the minimum A's also gets less and less smooth, as a increases. This 
phenomenon, together with the loss of hyperbolicity*, dramatically slows down 
the convergence of cycle expansions, as we discuss in the next section. Figure|31 
shows the trend of the eigenvalues in the 2B~ map when a = 0.1. As we can 
see, it is very difhcult to make a good fit for the trend of the minimum |A|'s 
in this case, because all the coupled T2 (and T3) maps have got almost stable 
cycles for such a's, including some at high period. It is interesting to compare 
this behaviour with what happens with intermittent ID maps ^22, . 

The other interesting issue we analyse is whether and why the number of 
cycles of a coupled map for a fixed period varies as a function of the coupling 
parameter a. Numerical simulation shows that few periodic orbits of the un- 
coupled map (O move off the fundamental domain when the system is coupled 
(a > 0). Moreover, the number of rip > 4 cycles continues to decrease as a 
grows, because some periodic solutions of the set (|5I6II become complex or de- 
generate for certain values of a, so that most of the missing orbits have just 

^henceforth we call A eigenvalue, but we always mean its real part. 

•^The same happens with the T3 coupled maps where the maximum [A] is 9^"p when a = 
and Q"?" or lower otherwise. 

*for a hyperbolic system |A| > Ae^"p for some A,b> 



|Ai|>l 
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Figure 2: Trend of maximum and minimum eigenvalues of cycles up to period 10 
of the 2A~^ map. At the top a = (uncoupled case), in the middle a = 0.0001, 
at the bottom a = 0.1. 
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Figure 3: maximum and minimum eigenvalues versus orbit period of the 2B 
map for a — 0.1. 



model 


V 


C 


2A+ 


0.50 


0.50 


2A- 


0.50 


0.53 


2B+ 


0.50 


0.50 


2B- 


0.50 


0.54 


3yl± 


0.50 


0.61 


3B± 


0.50 


0.63 



Table 1: Parameters -q and C of the fitting curve (O with a e [10 ^,10 
uncertainty of about 2 on the last digit 



disappeared into the complex plane. Data representing the number of periodic 
points H{a, Up) with a S [10^^, 10^'^] have been fitted with the curve 

H{a, Hp) = N^'"" exp(-2npCa'') (8) 

where C and rj are parameters that depend on the model employed, as shown in 
tabled As we can see, rj — 0.5 in all the models, which is theoretically plausible 
since the dynamics and bifurcations are controlled by the critical point of the 
Tchebyscheff map, which is quadratic. That implies that perturbations of order 
a have effects of order y/a, hence the exponent of 0.5. H{a, Up) is related to the 
topological entropy h{a) of the system by 

h^a)= lim i^feM 

rip^oo Up 

SO that 

h[a) = 2[lniV - Ca''] 

which decreases with a power law as a increases. Figure0]shows good agreement 
of the h{a) fitting function with the data for very weak coupling (a < 10^'^) 
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Figure 4: ln[21niV — h{a)] versus a in the 2^4+ model. 



while a slight deviation for larger a. This deviation is not unexpected, since 
this form for the topological entropy is only an approximate expression in the 
small coupling limit. Of course a detailed investigation could reveal that the 
topological entropy is not smooth in fine detail, however this form is remarkably 
accurate for practical purposes. 



3 cycle expansions and computation of averages 
3.1 cycle expansions 

Given a quantity v{x), we want to compute its average on the attractor of 
the dynamics (|5I6|) . such as the one depicted in figure [3 Assuming that such 
attractor is ergodic for any value of the coupling parameter we consider, direct 
simulation is the quickest method to use for this purpose, but its precision is 
limited by statistics. That is one reason why we use cycle expansions, which 
are potentially far more precise. According to periodic orbit theory [T^, the 
average of v can be found by computing the dynamical zeta function 

l/C = l[(^-tp) ,tp^^j^ (9) 
p \ p\ 

where 

n-l 
i=0 
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Figure 5: Scatter plot of {x,y) showing the attractor of the 2i?+ model at 
a = 0.1. 

/ is the map in question, /9 is a parameter, Ap the eigenvalue of the cycle p and n 
is the orbit period. The uncertain digit in the average is obtained by comparing 
results from the computation of lO up to the two largest periods (for example 
10 and 11). It is possible to check the rate of convergence of the zeta function 
by setting /3 = and z — 1 and computing 1/^(0,0), which should approach 
zero. This way we can have an idea of how effective our cycle expansion is and 
whether longer cycles are needed. Figure El shows the convergence of 1/^(0,0) 
versus the truncation period of the expansion. As we can see, the convergence is 
monotonic and exponential when the coupling parameter is relatively small (say 
a < 10^^), while it becomes nonmonotonic and increasingly slower for larger 
values of a. 

The accuracy of 1/C is particularly affected towards a — 0.1, where longer 
cycles or stability ordering may be needed. Ordering by stability means to 
expand the zeta function keeping all terms whose period and stability are less 
than a given cutoff (|A| < Amax)- This scheme is usually preferred to length 
ordering in case of intermittency or bad symbolic dynamics grammar |19| . but 
it hasn't been attempted here since it would need much a higher period of 
truncation for the cycle expansion and, as a consequence, impractical CPU 
times. A full resolution would require understanding which orbits are close to 
stability, which is beyond the scope of this paper. 

Anyway the behaviour shown in figure El traces back to the trend of the 
minimum eigenvalues of the cycles (sec. 2.3), since we can see from equation jnj 
that the smaller |A|'s are more relevant in the computation of the zeta function. 
That means that if their trend is exponential with respect to increasing a's. 
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Figure 6: Convergence of 1/C(0, 0) (logarithmic scale) versus maximum cy- 
cle length of the expansion, from the top left to the bottom right the maps 
3^1+ , 2_B+, 3i?+, shown for different values of the coupling a. 



the tp's are smaller and smaller as 1/C(0, 0) approaches zero. Otherwise the 
convergence gets more complicated. The expansion for the T3 maps has been 
truncated at period 7, yet the precision obtained is comparable with what we 
have got in the T2 maps. That is because 9*" « 4^^, meaning a similar number 
of periodic points and hence more or less the same number of terms in the two 
expansions. 



3.2 Computation of averages 

We now want to use cycle expansions to evaluate the average of an observable 
over the two-site lattice H5I6|I and compare the outcomes with what we get with 
direct simulation on the same lattice. According to the theory of Beck ^5], an 
interaction potential 

W{a) = (11) 

and self-interacting potentials 

y(^=3)(a) = -$4 + ^$2 (13) 
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Figure 7: Average of the self-interacting potentials V'^^^ (a) for the 2A+ and 2i?+ 
maps (left-hand side) and F'^^ (a) for the 3A+ and 35+ maps (right-hand side) 
versus coupling parameter a , using the expansion Q and direct simulation. 



can be defined for the chaotic strings described by equation Zeros of Hll|l 
and minima of Ij 121 13(1 may represent parameters of the Standard Model of 
particle physics such as coupling constants and masses of elementary particles. 
Since both the potentials defined above depend on a in a nontrivial way, we can 
only find zeros and minima numerically by averaging both V{a) and W{a) over 
the lattice (0J). Beck has obtained his results by running a long trajectory on a 
lattice with many sites, we are going to see what happens instead by computing 
the zeta function @ on the two-site lattice H5I6|I . Figure [7| shows our results 
(achieved by computing ^ up to two different orbit periods) compared with 
direct simulation of the potentials V'-^hA similar behaviour has been observed 
with the interaction potential W{a). Five main features emerge. 

1. In most points cycle expansion coincides with direct simulation. 

2. Cycle expansion is far offset on some points, which can be explained with 
the presence of almost stable cycles (| A| 1) of high periods whose terms 
in eq. Q have too much weight in the expansion. 

3. There are points where both cycle expansion and direct simulation are 
offset. That is probably caused by one or more stable cycles or chaotic 
attractors, whose basins of attraction alter the phase space making it no 
longer ergodic. 
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Figure 8: Average of the self-interacting potential y^^^ (a) in the 25+ model 
obtained with cycle expansion and direct simulation (with different starting 
values from those in figure EJ) versus coupling parameter a. 



4. The precision of the expansion visibly decreases as a approaches 0.1, as we 
expect from the observations made in the previous sections on the trend 
of the minimum eigenvalues and the behaviour of 1/C(0, 0). 

5. Figure IHl shows an offset at some point between cycle expansion and direct 
simulation in the model 2i?+ performed with different starting values from 
those used in figure [3 That proves a lack of ergodicity in some regions 
of the interval [0, 0.1] of a and indicates that the unstable periodic orbits 
are mostly unaffected by the presence of a stable orbit (or other ergodic 
component such as an elliptic region). Hence the cycles form a transient 
chaotic set more or less covering the attractor (Fig. O for similar values of 
the parameter. In this sense, the periodic orbit expansion results are more 
stable than direct simulation with respect to variation of the parameter, 
and may be preferred for this reason. 



4 Conclusions and outlook 

Periodic orbits have provided a unique window with which to view properties 
of weakly coupled CML. The periodic orbits in the uncoupled limit provide a 
reliable means of finding periodic orbits of the weakly coupled system. The 
symmetry of the coupled system has been exploited to improve convergence of 
the cycle expansion. Cycle expansion techniques work as expected, well near 
the uncoupled, more chaotic limit, but with slower convergence at stronger cou- 
pling. The computation of averages using cycle expansions leads to important 
information at a glance concerning the number and nature of nearly or fully 
stable states coexisting with the chaotic web of unstable periodic orbits. 

The observed power law dependence of the topological entropy with param- 
eter requires a more detailed investigation of the density and nature of bifurca- 
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tions in this system to explain the value of the power 77 = 1/2, the extent of its 
apparent universality and corrections for finite values of the parameter. For a 
more rigorous approach to the topological entropy of a different two site model, 
see 

This work has clarified a number of issues and challenges regarding a more 
general periodic orbit approach to CML. For stronger coupling (hence weaker 
chaos) stability ordering of cycle expansions is likely to be more effective 
than the length ordering employed here; this requires an understanding of the 
symbolic dynamics of the most stable cycles, which may be obtained in the study 
of the bifurcations mentioned above. Larger lattices are of course an important 
aspect, however the rapid decay of spatial correlations for weakly coupled CML 
make this less important than it might seem. The fact, as observed here, that 
small windows of stability at weak coupling affect the cycle expansions at only a 
few parameter values, gives great confidence since their presence is expected to 
be even less evident in lattices of larger size, in which typical orbits have many 
unstable directions. 
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